Measurement of the reaction ^^0(a,n)^°Ne and its impact on the s process in massive 

stars 
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Background: The ratio between the rates of the reactions ^^0(a, n)^''Ne and ^^0(a,7)^^Ne determines whether 
^''O is an efficient neutron poison for the s process in massive stars, or if most of the neutrons captured by ^®0(n, 7) 
are recycled into the stellar environment. This ratio is of particular relevance to constrain the s process yields of 
fast rotating massive stars at low metallicity. 

Purpose: Recent results on the (a, 7) channel have made it necessary to measure the (a, n) reaction more precisely 
and investigate the effect of the new data on s process nucleosynthesis in massive stars. 

Method: The ^^0(a, n(o+i)) reaction has been measured with a moderating neutron detector. In addition, the 
(a, ni) channel has been measured independently by observation of the characteristic 1633 keV 7-transition in ^°Ne. 
The reaction cross section was determined with a simultaneous i?-matrix fit to both channels. {ct,n) and (0,7) 
resonance strengths of states lying below the covered energy range were estimated using their known properties 
from the literature. 

Results: The reaction channels ^^0(a, no)^''Ne and ^^0(Q,ni7)^''Ne were measured in the energy range Eq, = 
800 keV to 2300 keV. A new ^^0(Q,n) reaction rate was deduced for the temperature range 0.1 GK to 10 GK. At 
typical He burning temperatures, the combination of the new (a, n) rate with a previously measured (a, 7) rate 
gives approximately the same ratio as current compilations. The influence on the nucleosynthesis of the s process 
in massive stars at low metallicity is discussed. 

Conclusions: It was found that in He burning conditions the (0,7) channel is strong enough to compete with the 
neutron channel. This leads to a less eflicient neutron recycling compared to a previous suggestion of a very weak 
(q, 7) channel. S process calculations using our rates confirm that massive rotating stars do play a significant role 
in the production of elements up to Sr, but they strongly reduce the s process contribution to heavier elements. 

PACS numbers: 26.20.Kn, 24.30.-v, 25.55.-e 



I. INTRODUCTION 

Most of the elements in the mass range 60 < A < 90 
that we observe today in the Solar system are produced 
by neutron capture on iron seed nuclei, mainly during the 
convective core helium and convective shell carbon phases 
in massive stars (weak s process) [1]. The efficiency of 
the weak s process depends on the network of nuclear 
reactions used for stellar calculations, above all on the 
main neutron source ^^Ne(a,n)^^Mg and on its main 
competing reaction ^^Ne(a,7)^^Mg, and on the set of 
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neutron capture cross sections used in the energy range 
of He and C burning (0.25 < T < 1.5 GK) [2-4]. 

Since the rate of neutron captures is slow compared 
to the decay rate of unstable reaction products, in the 
s process the neutron capture path follows the valley of 
stability. The final s process abundances depend on the 
total amount of neutrons available integrated over time (or 
neutron exposure) [5] , the history of the neutron density 
in the stellar regions where the s process occurs, and on 
the initial stellar metallicity. Light isotopes, depending on 
their abundance and neutron capture cross sections, can 
capture a large amount of free neutrons, thereby acting 
as neutron poisons in the burning environment. 

The main neutron source for the s process in massive 
stars is the ^^Ne(a, n)^^Mg reaction. It is activated mostly 
at the end of the convective core He burning [6] and in 
the convective carbon shell, on ashes of previous He core 
burning [7]. The abundance of ^^Ne is given by the initial 
CNO abundances, where "'^^N produced via H burning in 
the CNO cycle forms ^^Ne by capturing two a-particles 



in the initial He burning phases ([4, 8] and references 
therein). The s process in massive stars is a secondary 
process because of its dependence on the initial metallicity, 
the ^^Ne abundance, on the iron seeds and on the effect 
of light neutron poisons. 

Because the efficiency and yields of the s process de- 
crease linearly with the initial metal content for massive 
stars with lower mctallicities [9-11], the s process contri- 
bution to galactic chemical evolution from massive stars 
becomes marginal at low metallicity. However, recent 
observations (e.g., [12-15]) show a puzzling enhancement 
of light s process elements in very low metallicity stars. 
Recent theoretical studies have shown that fast rotating 
massive stars may potentially have s process yields orders 
of magnitude higher than non-rotating stars [16, 17]. In 
this case primary ^^Ne is produced in the convective He 
core, independently from the initial metallicity of the star. 
As a result, fast rotator yields could be relevant for the 
chemical evolution of heavy s process elements. 

Baraffe et al. [11] discussed the potential relevance 
of the ^^0(q;,7) and ^^0(a,n) rates for the s process 
in massive stars at low metallicity. Specifically, the rel- 
ative (a,n)/(a,7) ratio constrains the neutron poison 
efficiency of the ^^0(n,7)^^0 reaction. ^^O is the most 
abundant species in regions processed by advanced He 
burning via the ^^C(a, 7)"'^^0 reaction. Because of its high 
abundance -^^O is the strongest neutron absorber, despite 
the low ^^0(n,7)"'^^0 cross section [4]. The present large 
uncertainties of the ^'^O+a rates (together with other 
relevant rates, e. g., ^^Ne+a [16]) affects the theoretical 
s process predictions of fast rotating massive stars at low 
metallicity [18]. 

The neutron source ^^Ne(a, n)^^Mg is activated in mas- 
sive stars at stellar temperatures of typically T ^ 0.25 - 
0.3 and T ^ 1 GK, during convective core He burning 
and convective shell C burning, respectively. For the re- 
action ^^0+a, center-of-mass energies of around 0.5 MeV 
and 1.1 MeV are therefore most important. The reac- 
tion ^^0(a, n)^'^Ne has previously been measured in the 
energy range between 0.6 and 12.5 MeV [19-21]. The 
lower energy neutron data are only available in form of 
a Ph.D. thesis [19]. Its reliability is uncertain due to a 
strong background contribution from the -'^^0(a,n)^'^Ne 
channel, as well as the inability of the experimental setup 
to discriminate between the uq and ni channels, leading 
to ambiguities in the detector efficiency determination. 

To overcome these problems we remeasured the total 
yield of the reaction ^'^0(a, ntotai)^°Ne observed using 
a high efficiency neutron detector and targets with a 
very low (0.4%) ^*0 content. The experimental setup 
and the measurement are described in sections II B and 
IIIB, respectively. In an additional experiment the 
^^0(a,ni7)^''Ne channel was measured by observation 
of the E = 1633 keV 7-transition to the ground state of 
^°Ne (sections HA and III A). This allowed a subtraction 
of the ni contribution to the total yield and a better 
determination of the effective neutron detector efficiency. 
In section IV B we present a simultaneous i?- matrix fit to 



both neutron data sets. 

Until recently the only available rates for the com- 
peting ^^0(a,7)^^Ne channel were an estimate from the 
Caughlan- Fowler compilation [22] (CF88 hereafter) and 
a rate on the basis of microscopic cluster model calcula- 
tions [23]. These two rates differ by a factor of ~ 1000, 
causing a large uncertainty in the s process yields from 
fast rotating massive stars [18]. A first measurement of 
the i'^0(a,7)2iNe reaction [24] supports the CF88 esti- 
mate. This recent experimental data on the (a, 7) channel 
and the (a, n) data presented here were combined with 
low-energy extrapolations and used as input for a stel- 
lar network calculation. The contributions to both the 
(a, 7) and (a,n) channels from states below the energy 
range covered in this experiment are estimated in sec- 
tion IV C. The new reaction rates and a discussion of the 
astrophysical implications can be found in section IV D. 



II. EXPERIMENTAL SETUP 

The experiment consisted of two independent measure- 
ments, a direct measurement of the reaction neutrons and 
the detection of the characteristic 1633 keV 7-ray from 
the "^^0(a, ni7)^''Ne channel. The specific setups for each 
part are described separately in the following subsections. 
Here we describe the parts of the setup that were common 
to both measurements. 

The a-beam was provided by the 4 MV KN acceler- 
ator at the University of Notre Dame Nuclear Science 
Laboratory. Energy calibration and resolution (1.1 keV) 
were determined using the well-known resonances Ej, = 



991.86 ± 0.03 keV and E„ 



1317.14 ± 0.07 keV in 



^^Al(p,7)^®Si [25]. The beam energy was reproducible 
within ±2 keV between different hysteresis cycles of the 
analyzing magnet during the course of the experiment. 

Targets were prepared by anodization of 0.3 mm thick 
tantalum backings using H2O enriched to 90.1% in I'^'O**. 
This process is known to produce homogeneous films of 
Ta205 [26, 27]. The film thickness can be controlled in 
a reproducible way through regulation of the maximal 
anodization voltage. The target thickness was chosen to 
be about 9 keV of energy loss for an a beam of 1000 keV 
and a target orientation of 90° with respect to the beam 
direction. 

To reduce carbon deposition a liquid nitrogen cooled 
copper tube (cold finger) was mounted in front of the 
target. A bias of -400 V was applied to the cold finger 
for suppression of secondary electrons. The beam was 
scanned with magnetic steerers to produce a beam spot 
size of 1.4 cm x 1.6 cm on the target. In both parts of the 
experiment the target chamber was electrically isolated 
for charge collection and the targets were directly water 
cooled using deionized water. 



Purchased from Isotech, Miamisburg, OH. The ^*0 and ^^O 
contents of the water were 0.4% and 9.5 %, respectively. 
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FIG. 1. (Color online) Efficiency of the germanium detector 
as a function of photon energy. The solid line represents a fit 
to the data points. 



A. ^^0(a,ni7)^"Ne setup 

For this part of the measurement a Ge detector with a 
relative efficiency of 20% was placed at a distance of 5.4 
cm from the target at an angle of 45° with respect to the 
beam direction. In order to minimize radiation damage 
to the detector a 2.5 cm polyethylene disk was attached 
to the detector front cap. This reduced significantly the 
neutron flux in the detector. The transition from the 
first excited state to the ground state of ^°Ne emits a 
7 ray with an energy of E^ = 1633.7 keV [28]. Abso- 
lute efficiencies were established with calibrated ^'^^Cs, 
^°Co and ^^'^Ba sources and augmented using relative 
efficiency data from ^^Co. There was enough absorbing 
material between source and detector (0.25 mm tantalum 
backing, 1 mm brass target holder and the 2.5 cm thick 
polyethylene disk) to attenuate the low-energy photons 
emitted by the Ba source (E^ = 81 keV and E-y — 80 keV) 
by approximately 98%, resulting in negligible summing 
effects in the Ba measurements. Summing corrections for 
the ^°Co measurements amounted to less than 2%. The 
resulting efficiency curve is shown in Fig. 1. 



B. 



'0(a,ntota!)^°Ne setup 



This part of the experiment was dedicated to the mea- 
surement of the ^^0(a,ntotai)^*'Ne reaction by direct de- 
tection of the reaction neutrons. The detector consisted 
of 20 "^He counters that were embedded in an arrangement 
of two concentric rings with diameter 6 cm and 11 cm, 
respectively into a 30 cm x 30 cm x 33 cm polyethylene 
moderator. The very high cross section of the reaction 
3He(n, p)3h (Q = 764 keV, a = 5330 barn) for thermal 
neutron capture and its low sensitivity to photons make 
it an excellent neutron detector, although due to the ther- 
malization in the polyethylene moderator the information 
on the neutron energy is mostly lost. 



FIG. 2. (Color online) Geometry of the neutron detector used 
in the Geant4 simulation. Shown are the outlines of the outer 
borated polyethylene shield, the moderator with the 20 '^He 
counters and the beam line with the target holder in its center. 
Depicted in black is the simulated track of a neutron. 
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FIG. 3. (Color online) Efficiency of the neutron detector 
as a function of energy. Shown are the results of activation 
measurements and the Geant4 simulation. Above a neutron 
energy of E„ = 650 keV the cross section of the ^^V(p,n)^^Cr 
reaction is high enough that dead time effects prohibited 
further measurements. The MCNP data are not shown here 
but follow the shape of the Geant4 results. 



Additional shielding from background neutrons was 
provided by a 5 cm thick layer of borated polyethylene 
on the outside of the moderator. At the beginning of 
each experimental run the detector was centered around 
the target chamber by monitoring the neutron yield as a 
function of distance of the detector to the target holder. 
The position with the maximum yield was then used as 
the default position for all further measurements. The 
target was mounted at an angle of 90° with respect to 
the beam direction. 

The absolute detector efficiency was determined using 
the reaction ^^V(p,n)^^Cr [29-31]. This reaction with a 
neutron threshold of 1565 keV is known to produce mono- 



energetic neutrons that are approximately isotropic in the 
center-of-mass system. The thickness of the vanadium 
targets was about 30 [igjcm} ^ corresponding to a proton 
energy loss between 2.5 and 3 keV in the covered energy 
range. Although a large number of narrow resonances 
has been observed in this reaction no pronounced angular 
distributions have been reported in the literature [32, 33]. 
Since our data show a smooth behaviour, in agreement 
with similar calibration measurements using ^^V(p, nf^Cr 
[34, 35], we believe it to be unaffected by any possible 
angular anisotropics. 

The reaction product ^^Cr decays to the first excited 
state in ^^V by electron capture and subsequent emission 
of a characteristic E^ = 320 keV photon with a half-life of 
Ti/2 = 27.7 d and a branching probability of 9.9 % [36]. 
The activity and therefore the number of produced ^^Cr 
nuclei was determined by counting the emitted 7-rays 
using a Pb-Cu shielded counting station at Notre Dame. 
In this way the neutron detection efficiency of the detector 
was established with high accuracy without having to 
rely on cross section data. The experimental uncertainty 
is dominated by the efficiency determination of the Ge 
detector (2%) introduced by the uncertainty in the source 
activity and errors due to the activation procedure. To 
obtain sufficient ^^V activity it was necessary to run 
with large beam currents for which the count rate in the 
neutron detector was too high. Therefore, the number 
of generated neutrons per accumulated charge had to be 
determined during several short runs with reduced beam 
current. This number was then used to scale the higher- 
current activation runs. We estimate the error resulting 
from this procedure to be ±5%. The resulting total error 
of the activation measurement is 5.4%. 

The ^^V activations were carried out at proton energies 
between 1600 keV to 2200 keV resulting in calibration 
points for neutron energies between 50 keV and 650 keV. 

In addition to measurements Geant4 [37] and MCNP 
[38] simulations of the setup were used to determine the en- 
ergy dependence of the efhciency towards higher neutron 
energies that were inaccessible to direct measurements 
[39]. The geometry that was used for the simulations 
is shown in Fig. 2 and the resulting efficiency can be 
seen together with the measurements in Fig. 3. A more 
detailed description of the detector and the efficiency de- 
termination can be found in Ref. [40]. The simulations 
assumed an isotropic angular distribution of the emitted 
neutrons. Due to the large dimensions and almost com- 
plete 47r coverage of the detector the effects of possible 
anisotropics are strongly reduced in our measurements. 



C. Target thickness 

The observed yield F of a nuclear reaction as a function 
of energy can be described by [41]: 
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FIG. 4. Yield (reactions per projectile) of the Eq = 1002 keV 
resonance in ^'^0(Qf, 7)^'^Ne that was used for the determination 
of the target thickness. The line is a fit to the data points. 



where Eq and A_E stand for the energy of the incident 
particle and the target thickness in terms of energy loss; 
CT and e are the energy-dependent reaction cross section 
and the effective stopping power of the target material. 
Due to the statistical nature of energy loss a second 
integration over a distribution /(i?, i?') that describes 
the probability for a projectile having energy E' at depth 
E in the target can be included in the yield calculation. 
Under the assumption of Bohr straggling, f{E,E') can 
be approximated as a Gaussian with an energy-dependent 
width [42] . The target thickness AE and the range of the 
yield integration is related to the number density of active 
(^^O) target atoms n through AE — ne [43]. Therefore, 
the target thickness needs to be known for extracting the 
cross section from yield data. 

The relationship between the anodization voltage 
{AU = 5 V was used for the production of our tar- 
gets) and the thickness of the resulting Ta205 layer is 
t = 1.9225(y + 1.4) nm [27]. The expected thickness 
of 12.3 nm translates into a nominal area density of 
n = 6.9 X iQie oxygen atoms ^ rj.^ thickucsses of various 

targets were determined experimentally from the yields 
of isolated resonances in the reactions ^"^0(0, 7)^"'^ Ne (Eq, 
= 1002 keV), "O(a,n)20Ne (E^ = 1247 and 1841 keV) 
and i^O(a,ni)20Ne (E„ = 1841 and 2020 keV). 

Fig. 4 shows the measured reaction yield of the 
Ea = 1002 keV resonance in ^'^0{a,-/)^^Ne. The line 
through the data points is a fit of the yield based on 
equation 1, where the single-level Breit-Wigner formula 
"' o^ (•j;-j;°-)2+r^/4 ^^ used to describe the energy depen- 
dence of the cross section. The average ^^O area density 



Y = 



a{E' 



k,-AE Jo ^i^' 



j{fiE,E')dE'dE, (1) 



resulting from the above mentioned 5 resonance scans is 
n = (7.5 ±0.4) X 1016M5S5, Sjnce a number of targets was 
used over the course of the experiment a slight variation 
in the individual target thicknesses is to be expected. We 
varied the input thicknesses for the i?-matrix calculation 
and found the minimum x^ value of the fits at an area 
density of n = 7.8 x iQie atoms ^ r^j^j^ demonstrates that 
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FIG. 5. Yield (reactions per projectile) of the ^^0{a, ni7)'^°Ne 
reaction. Also shown is the threshold of the rii channel at 
Ea ~ 1293 keV. The arrows denote upper limits. 



the fluctuation of the target thicknesses is within the 
quoted error. 



III. EXPERIMENTAL RESULTS 
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FIG. 6. (Color online) Yield of the ^^0(a,ntota;)^°Ne reaction. 
Also shown in the lower curve in red is the neutron yield 
obtained with a blank tantalum target. The line through 
the blank yield is to guide the eye only. Clearly visible in 
both curves are the strong 1054 and 1336 keV resonances 



^C(c 



)^^0. The error bars show the statistical error 



only. Above E^ = 1293 keV both the no and ni channels 
are open and contribute to the yield. The dashed line is the 
^^C(a,n)'^^0 cross section from Ref. [44]. It has been scaled 
for comparison with the yield data. 



A. "0(Q,ni7)2°Ne 

An excitation curve of the ^^0(a,ni7)^°Ne reaction 
from the m threshold at £„ = 1294 keV to 2300 keV was 
measured in steps of 5 keV or less. For this measurement 
the E..y = 1633.7 keV transition from the first excited state 
to the ground state in the ^°Ne nucleus was observed. 
Angular correlation effects between the direction defined 
by the incident beam and the detected photons were 
calculated to be less than 8 % in a worst-case scenario. 
Since the emitted neutron is unobserved in this part of 
the experiment and the spin of the intermediary state is 
unknown angular distribution coefficients were calculated 
for various spin assumptions for the intermediary state 
in ^^Ne. The yield Y (number of reactions per projectile) 
was calculated from the intensity / in the 1633 keV peak 
by: 



Y 



Qdtr\ 



(2) 



Qdt and r\ represent the dead time corrected number of 
projectiles and the detector efficiency at 1633 keV. The 
resulting excitation curve is shown in Fig. 5. 



B. i'0(Q,n,,i„i)"'Ne 

Fig. 6 shows the ^'^ 0{a,ntotaif'^^^ reaction yield mea- 
sured with the neutron detector described in Sec. II B. 



The yield was not corrected for the detector efficiency. At 
a-energies above 1.3 MeV population of the first excited 
state in ^°Ne is energetically possible and two neutron 
groups contribute to the yield. The 1.05 MeV difference 
in energy between neutrons from the two groups, and the 
inability to discriminate between the two channels with 
the moderating detector, results in an uncertainty in the 
detector efficiency. This problem will be discussed in the 
following sections. 

Information on tlie initial energy of the neutron is 
lost during the moderation process in a detector of the 
type used in this part of the experiment. Therefore, 
background neutrons are in principle indistinguishable 
from the reaction neutrons of interest. To avoid the 
assignment of spurious resonances during the analysis of 
the data one must carefully investigate possible beam- 
induced neutron producing reactions. The most important 
background reaction is ^■^C(a,n)^^0 due to its liigh cross 
section and inevitable carbon buildup on beam collimators 
and slits. It has two strong resonances at Eq = 1054 keV 
and 1336 keV [44] that lie within our energy region of 
interest. To investigate the beam-induced background a 
yield curve was measured with a blank Ta target, covering 
the whole energy range of the experiment. The result 
of this is shown in the lower curve in Fig. 6. The two 
^•^C resonances are clearly visible in both the ^^O and 
the blank data. The sensitivity limit has been reached 
towards the lowest energy points of our measurement, and 
the beam-induced background dominates the yield. For 
comparison the ^'^C(a, n)^^0 cross section is also shown in 



Fig. 6 (scaled by a constant factor). It can be clearly seen 
that ^"^C is indeed the source of background neutrons over 
the whole energy range. The typical amount of carbon can 
be estimated using the thin-target yield Y = an, where Y 
is the backing yield and a the ^^C(a,nY^O cross section. 
The resulting "c density is n « 10" s^. The lower 
energy ^^C resonance was also used as an additional beam 
energy calibration point for our experiment. 



IV. ANALYSIS 
A. Separation of the (a, no) channel 

Because of the occurrence of two neutron groups 
with different energies above the ni threshold at Eq, — 
1294 keV, the neutron detector efficiency can only be 
determined accurately if the branching of the Uq and ni 
reaction channels is known. Since we independently mea- 
sured the ^''0(a,ni7)^''Ne reaction by detection of the 
characteristic E.y — 1633.7 keV transition, its contribu- 
tion to the ^^0(a, ntota/)^°Ne reaction could be subtracted 
from the data. As the target thicknesses used in the two 
experiments were different due to the different angle of 
the target with respect to the beam axis (90° and 45°), 
we determined the cross section of the ^^0(a, ni7)^°Ne 
reaction channel with a preliminary i?-matrix fit. The 
simultaneous i?-matrix analysis of both channels and the 
computer code used for the calculations is described in 
the following Section IV B. 

Fig. 7 shows the results of the preliminary i?-matrix fit 
to the (a,ni7) data (with degree, of freedom = ^-O)- It was 
necessary to include a background pole that destructively 
interferes with the £„ = 2150 keV resonance. As there is 
little information on the J'^ values of the included states 
available in the literature, random spins corresponding to 
low angular momenta were assigned. The only restriction 
was that states that clearly interfered with each other had 
to be assigned the same spin-parity value. 

The R matrix cross section of the {a,ni) channel cr„j 
was then converted into an "experimental" yield contribut- 
ing to the utotai measurement (at a energy Eq) ^^(-E'o) 
by integration of the cross section over the target thick- 
ness AE and taking into account the efficiency of the 
neutron detector ri{En^) and the stopping power of the 
target material e(£^o)- 90 degree neutron energies were 
used in this and all later calculations. 
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FIG. 7. (Color online) Preliminary 7?-matrix fit to the 
^^0(Q,ni7)^''Ne data. The line is the fit to the data points. 
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FIG. 8. Yield of the "0(a,no)^''Ne reaction channel. For 
this plot the integrated ^^0(Q,ni7)^°Ne cross section was 
scaled with the neutron detector efficiency and subtracted from 
the ^'^0{a,ntotaiY'^"^^ yield. The yield in the energy regions 
around Eq = 1050 keV and E^ = 1350 keV is dominated by 
the two low-energy resonances in ^''C(Q,n)^^0. These data 
points and the data in a region around 1200 keV where a 
^*0(a, n) resonance contributes to the yield were removed for 
the final analysis. The arrows denote upper limits. 



cancels out by the symmetric setup about 90 degrees and 
by the very large angular coverage of the detector. 



YnAEo) 



E„ 



En-AE 



dE^rjiEr.,] 



(3) 



This contribution was subtracted from the measured 
total neutron yield Ytotai and the result scaled with the 
efficiency for the neutron energy of the uq channel to give 
the separated no yield Y„g shown in Fig. 8: Yno{E) — 
v{Eno){YtotaiiE) - Y„^{E)). Any significant effect of the 
neutron energy distribution on the detector efficiency 



B. U-matrix calculations 

The i?-matrix analysis was performed in the framework 
of a multi-level, multi-channel approach based on the 
formalism outlined for the i?-matrix code AZURE [45]. 
The multi-channel capabilities of the code enabled us to 
simultaneously fit both measured reaction channels. In 
order to directly fit the i?-matrix yield to the measured 
yield data a target integration routine was added to the 
program. The neutron yield from the blank Ta target 
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FIG. 9. (Color online) Experimental and calculated yield 
of the ^'^0(Q,no)^''Ne and ^'^0(Q:,ni7)^'^Ne reaction channels. 
The lines through the data points are the result of an ij-matrix 
fit to both channels. The ni channel has been divided by a 
factor of 10 to better separate it from the no plot. The arrows 
denote upper limits. 



FIG. 10. (Color online) S-factors of the "0{a,nof°Ne and 
^^0(a,ni7)^°Ne reaction channels. Shown as the points is the 
deconvolved yield data (see text) for both channels and the 
respective R matrix s-factors are represented by lines. The 
s-factors of the (a, ni'y) channel were multiplied by a factor 
of 0.1 to enhance the visual separation of the curves. 



(see Fig. 6 was subtracted from the ntotui data before the 
final fit was conducted. The iZ-matrix parameterization 
following Brune [46] was used in the program. The channel 
radii were set to 5 fm and each target integration was 
divided into 25 sub-points. It should be emphasized 
that because of the ambiguity in the spin assignment the 
best-fit cross section can very likely be reproduced with 
different J'^ values and partial widths for the individual 
resonances. Therefore, the R matrix fit should be mainly 
considered as a deconvolution of the yield data and only 
the cross section and the resonance energies should be 
regarded as physically meaningful parameters. 

The measured yield, a parametrisation of the stopping 
power of the target e (tabulated e values were calculated 
with the computer code SRIM-2010 [47]) and the number 
density of ^^O atoms in the target (as determined in Sec. 
II C) were used as input for the ii-matrix calculation. 
For each yield data point (at energy Eq) the program 
calculated a cross section a, and used the integrated yield 
(Eq. 1) and the measured data for a least-squares fit with 
the Migrad algorithm from the Minuit2 minimization 
library [48] . It turned out that energy straggling of the 
projectile in the relatively thin target did not play a role 
and did not have to be included in the final calculation. 

The same background pole as for the ni data was 
included in the calculation. The 1.695 MeV and 1.84 MeV 
resonances are only observed in the ng channel. The 
measured yields and the results of the i?-matrix fit to 
both data sets is shown in Fig. 9. To make it easier to 
distinguish between the two channels the ui plot has been 
scaled down by a factor of 10. The reduced x^ values of the 
fit are 2.85 for the ni and 3.75 for the uq channel. The best- 
fit i?-matrix parameters from our calculation are shown in 
Table II. The relationship between the "observed" partial 
widths F and R matrix reduced widths 7 (as defined in 



Ref. [49]) is F,, 



2Pc7-e 



[46], where Pc and 



Sc are the penetration and the shift factor. Also listed 
are the energies of previously observed states from Ref. 
[25] . Some of these states have been seen as resonances 
in ^^0(a, n) [20] and ^°Ne(n,n) [50]. Known literature 
widths are listed for comparison. Values for the spin-parity 
were only available in two cases (at E^ — 8.065 MeV and 
E^ = 8.68 MeV). 

The calculated s-factors for both channels are shown 
in Fig. 10. Also shown in Fig. 10 are the experimental 
s-factors for both channels that have been obtained by 
an least-squares deconvolution method similar to the one 
described in Ref. [51]: in short, yield data was calculated 
by integration of the cross section over the target thickness. 
The starting point was the cross section obtained from 
the R matrix fit, which was then varied until the new 
integrated cross section agreed with the experimental 
yield. 

The uncertainty in the absolute scale of the s-factor is 
mostly influenced by the efficiency determination of the 
neutron and 7-ray detectors, and by the target thickness 
that is used in the integration process. i?-matrix fits for 
various thicknesses within the experimentally determined 
range (section II C) were performed. The effect on the 
calculated s-factor was ±10%. Following systematic un- 
certainties have to be attributed to the data: The error in 
the (a,ni) s-factor was calculated by quadratic addition 
of the ±10% from the target thickness, ±8% from our 
estimate of the maximal effect of the angular correlation 
and ±5% in the efficiency determination to be ±14%. The 
uncertainty in the lower-energy (a, no) data below the ni 
threshold (E^ < 1.3 MeV) is ±11% (10 % target thickness 
and 5.5 % detector efficiency). Finally, since the ng data 
above the ni threshold is the result of a combination of 



both 7 and neutron measurements the associated error is 
the quadratic addition of both uncertainties, or ±18%. 



C. LOW-ENERGY ESTIMATE 

The data presented here only extend down to E^^'' w 
800 keV. To cover the whole astrophysically relevant range 
it is necessary to estimate the contribution of lower-energy 
states. 13 states are known to exist between our lowest 
energy measurement and the a threshold in ^^Ne [25, 52]. 
Their possible contribution was evaluated as follows: At 
low energies charged-particle widths are generally much 
smaller than neutron and 7 widths (Ta <C ^n,-y)- Using 
an estimate of the a width and the branching between the 
7 and neutron channels r = p^ of a state the resonance 
strengths uijn,-/ can be evaluated from: 



^l{a,n) = ^ 



^l{an) = ^ 



Ljli^'yi^a 



T r 



r + r + r 



UjTn 



UjTa 



r„(r + i) 
r^(r-i + 1) 



r + 1 



(4) 



1 



(5) 



The required a widths can be calculated from alpha- 



particle widths Ta^sp and the a spectroscopic factor of 

the respective state: Fq, = C SFq,, sp- The computer code 

DWUCK4 [53] was used to evaluate the single-particle 

widths using a Woods-Saxon potential with radius rg = 

1.25 fm and diffuseness a = 0.6 fm. A spectroscopic factor 

of 0.01 was assumed for all states because no experimental 

information is available, 
r 
Information on the p^ ratios was inferred from the 

available literature: The relevant states have been seen 
in neutron scattering [50], ^^C(^^C, a)^^Ne [54] and 
i2C("C,a7)2iNe [55], (^He, p) [56] and ^^0{a,n-f)^^Ne 
[57]. A state at Ex = 7.96 MeV has been observed in 
^^0(^Li, np) [58]. Levels with a known width were as- 
sumed to have a strongly suppressed 7 channel and were 
assigned r = 10~^. The typical ratio between gamma and 
neutron widths of 10""^ was used if only the energy of the 
state is known. In the case of observed 7-deexcitation 
assignments between 0.1 and p^ = 2 were made based on 
the observed transition strengths. 

Our estimates and the resulting resonance strengths 
of the low-energy states used in our extrapolation are 
listed in Tab. I. To obtain an upper limit on their 
strengths, states with unknown spin and parities were 
assumed to be s-wave resonances. The lowest-lying state 
{Ex = 7.36 MeV) at only 12 keV above the a threshold 
cannot contribute to the reaction rate and was ignored in 
our extrapolation. 



TABLE L ^''O+a single-particle widths and estimated reso- 
nance strengths. E^ and 2J'^ from Refs. [25, 52] 


E^ [MeV] 


Bern [MeV] 


2.r 


1 


Ts, [eV] 


I7 
r„ 


^1-1 [^eV] 


a;7„ [^eV] 


7.422 


0.074 


(7+,ll-) 


2 


1.7 X IQ-'^'^' 


0.1 


2.1 X 10"^" 


2.1 X 10-2^ 


7.470 


0.122 


(1,3)- 


1 


1.5 X 10-22 


10"^ 


5.1 X 10-22 


5.1 X 10-^^ 


7.556 


0.208 







6.0 X 10"" 


10"^ 


6.0 X 10"" 


6.0 X 10-^° 


7.600 


0.252 







1.8 X 10"" 


0.1 


1.6 X 10"® 


1.6 X IQ-'^ 


7.627 


0.279 


3- 


1 


1.4 X 10"^'' 


10"^ 


9.5 X 10-^2 


9.5 X 10-^ 


7.653 


0.305 


7+ (5+) 





3.0 X 10"^ 


2 


4.0 X 10-5 


1.3 X 10-5 


7.750 


0.402 







2.2 X 10"^ 


10"^ 


2.2 X 10-5 


2.2 X 10-2 


7.815 


0.467 







5.4 X 10~^ 


10-^ 


5.4 X 10-'' 


0.54 


7.96 


0.612 


11" 


3 


1.2 X 10-"' 


2 


0.24 


0.08 


7.979 


0.631 


3- 


1 


8.3 X 10^3 


10-^ 


5.5 X 10-" 


55 


7.981 
8.008 


0.633 
0.66 


(7,11)+ 
1" 


2 
3 


1.3 X 10"^ 
3.1 X 10"'' 


1 
10-5 


8.9 
1.0 X 10-5 


8.9 
1.0 



The (a, n)) reaction rate (in units of cm'^mol ^s "') was 
calculated from the cross sections using the expression 

[59]: 



Na{'Jv)oi 



2.074 X 10 

T3/2 



10 



Ea{E)e 



-11.605B 



'^dE 



(6) 

The ^■^0(a,7)2^Ne rate was adopted from Ref. [24]. 
Both the (a, 7) and (a^n) rates were extrapolated to 
lower temperatures using the equation for isolated, narrow 
resonances: 



Na{(Jv) = 



2.643 X 10 

T3/2 



10 






(^7)^6 



-11.605£i/T 



(7) 



where resonance strengths uj^ were taken from Tab. I. 

For both equations the temperature T is given in GK 
and the cm. energies Ei and resonance strengths (^7)^ 
in MeV. The reaction rates are listed in Tab. III. The 
ratio of our recommended (a,?^) and (0,7) rates to our 
experimental rates is shown in the upper half of Fig. 11 
and 12 respectively. Shown in the lower half of the figures 
is a comparison of the recommended rates to literature 
rates. In the case of the {a,n) channel, the rate is shown 
in comparison to the more recent NACRE compilation 
which is based on unpublished experimental low energy 
data, and supersedes the previous CF88 rate. For the 
(01,7) channel our rate is compared to the CF88 rate 
because NACRE does not list a rate for this reaction. 
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FIG. 11. (Color online) Upper plot (a): Ratio of the recom- 
mended ^^0(a,n) rate to the experimental rate. Lower plot 
(b): Ratio of recommended ^'^©(a, n) rate to NACRE. 
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FIG. 12. (Color online) Upper plot (a): Ratio of the recom- 
mended ^^0(a,7) rate to the experimental rate. Lower plot 
(b): Ratio of the recommended ^'^0(a,7) rate to CF88. 



Plotted in Fig. 13 is the ratio of the NACRE {a,n) to 
the CF88 (01,7) rate, as well as the same ratio using the 
recommended rates from this work. At the temperatures 
of relevance for the s process, these two ratios are within 
a factor of 2. 



D. Astrophysical implications 

As discussed in the introduction, the 
^^0(a,n)/^^0(Q;,7) reaction rate ratio determines 
how efhciently ^^O acts as a neutron poison for the 
s process in massive stars. Uncertainty in this ratio 
predominantly arose from a lack of knowledge regarding 
the strength of the (0,7) channel. In previous estimations 
the ratio ranges from between « 10 (if using the CF88 
rates [22]) to « 10000 (if [23] is used). 

In the present work, both the (a, n) and the (a, 7) rates 
are weaker than those proposed by NACRE and CF88, 
respectively, at He burning temperatures. For instance, 
at T ~ 0.3 GK, the recommended (a, n) rate is about 1/4 
of NACRE, and the (a, 7) rate is about 1/4 of the CF88 
rate. 

Shown in Fig. 14 is the effect of our recommended 
-^^0(a,7) and ^^0(a,n) rates (red stars) on s process el- 
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FIG. 13. (Color online) The (a, n)/{a, 7) ratio as a function of 
temperature using the recommended rates, and compilations. 



emental abundances at the end of He core burning for 
fast rotating massive stars. This is compared to the abun- 
dances obtained using the ^^0(a,n) rate adopted from 
NACRE [60] and the ^^0(a,7) rate from CF88, modified 
by Descouvemont (blue squares). In both cases, the s 
process elemental abundances are plotted relative to the 
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FIG. 14. (Color online) Comparison of s process elemental abundances obtained at the end of convective He-core burning, 
relative to the initial abundance distribution. Red stars are the s process abundances using the ^^0(q, 7)^^Ne and ^^0{a, n)^'^Ne 
rates obtained in this work (where the ratio between the two rates is « 10 at T=0.3 GK). Blue squares are abundances obtained 
using the "0(a,7)^^Ne CF88 rate as modified by Descouvemont et al. [23], and the "0(a, n)^°Ne rate from NACRE (where 
the ratio between the rates is « 10000 at T=0.3 GK). See text for details. 



initial abundance. For comparison to Fig. 1 of Ref. [17], 
Fig. 14 is plotted as a function of atomic number. We also 
performed the calculations using the (a, 7) rate of CF88 
and the (a, n) rate of NACRE. The resulting elemental 
abundances are very close to the results obtained with 
the present rates, because both sets of rates have similar 
(a,n)/(a,7) ratios (see Fig. 13). 

Nucleosynthesis calculations were performed using the 
post-processing network code PPN [18]. The single-zone 
trajectory was extracted from a complete 25 solar mass 
star (e.g., Hirschi [61]), calculated using the Geneva stellar 
evolution code (GENEC) which is described in detail in 
Refs. [62] and [63], and assuming an initial metallicity 
of Z = 10~^. In order to take into account the primary 
•^^Ne production arising from stellar rotational mixing [16, 
17], the initial abundance of ^^Ne was set to be 1%, in 
accordance with Hirschi et al. [64]. 

Figure 14 demonstrates the sensitivity of the s process 
in low-metallicity stars to the ratio of the (a,n)/{a,j) 
channels. Our new measurement of the ratio confirms that 
the s process is not suppressed in fast-rotating massive 
stars, and that the main region receiving the additional 
feeding is between the Fe seeds up to Ba, peaking with 
maximum efficiency in the Sr region. Compared to the 
Descouvemont-modified CF88 ^^0(a,7) ratio however, 
we observe in general less s process efficiency because the 
modified (a, 7) rate is many orders of magnitude smaller 
than our measurement. As a result the neutron recycling 
is far more efficient. 



V. CONCLUSION 

We have presented new experimental results and 
i?-matrix calculations for the ^^0(a, n)^''Ne and the 
^^0(Q;,rii7)^°Ne reactions. A reaction rate for the {a,n) 
and the competing (a, 7) channel has been calculated 
and extended to lower temperatures by estimates on the 
contributions of resonances in lower-lying states. Despite 
having lower {a,Ti) and (0,7) rates than NACRE and 
CF88, the new recommended ratio has a similar impact 
on the s process elemental abundances. This is because 
the difference between the two ratios is less than a factor 
of 2 in the temperature range of relevance (see Fig. 13). 
Such minor changes would not be apparent on a log scale 
as used in Fig. 14. Using the CF88 (a, 7) and the NACRE 
(a, n) rates for a star with an average equatorial velocity of 
{veq) — 428 km s~^, Ref. [17] finds an s process elemental 
distribution which is similar to this work. Therefore, we 
confirm that massive rotating stars do play a significant 
role in the production of elements up to the Sr neutron 
magic peak. 

However significant uncertainties still arise from unob- 
served low energy resonances. A calculation to test the 
impact of the low energy resonances strengths shows that 
the nucleosynthesis is not sensitive to the absolute mag- 
nitude of the resonances, but rather to the relative ratio 
of the (a,7)/(a,n) reaction channels. Changing the ratio 
in the favor of (a,n) by a factor of 3 results in significant 



11 



enhancement of the A = 40 - 60 region. 

These results indicate that ^^O is a significant neutron 
poison in metal poor massive rotating stars. An extension 
of experimental data towards lower energies can further 
improve the precision of the reaction rates and seems 
feasible. 
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TABLE II: Best-fit _R-matrix parameters for the fits shown in 
figures 9 and 10. 1 and s are the relative angular momentum 
and the channel spin and F is the "observed" partial width as 
defined in Brune [46]. The "-I-/— " column shows the relative 
interference signs needed to reproduce our result. The states 
are sorted by ascending total angular momentum J. 

E^ [MeV] J'^ Channel 1 s r]eV] -h/- E^ (lit.) [25] 



{Continued.] 



E^ [MeV] J'' Channel 1 s T [eV] +/- E^ (lit.) [25] 



8.791 



8.189 



8.292 



8.438 



172^ 
3/2" 
3/2" 
3/2" 



;.47 3/2" 



.665 3/2" 



8.899 



9.203 



3/2" 
3/2" 



8.069 3/2+ 



.146 



.359 



8.839 



3/2+ 
3/2+ 
3/2+ 



no 
a 
no 
a 
no 
a 
ni 

no 
a 
ni 

no 
a 
no 



no 
a 
ni 

no 
a 
no 

a 
no 
a 
no 
a 
ni 



2 5/2 

1/2 

1 5/2 
1/2 
5/2 
1/2 
5/2 

1 3/2 
1 5/2 
1 1/2 
1 5/2 
1 3/2 
1 5/2 
1 1/2 
1 5/2 



627 

342 
88.4 xlO~^ 
3.32 xlO^ 
559 xlO-3 
7.54 xlO^ 
563 xlO"^ 
1.19 xlO^ 

140 
2.34 xlO^ 
363 xlO"^ 

320 

37.7 

781 

76.1 



1 1/2 62.5 xlO'^ 



531 
75.2 xlO^ 
1.83 xlO^ 

758 

85.6 
1.03 xlO^ 



5/2 

1/2 

5/2 

3/2 

5/2 

1/2 

5/2 46.2 xlO"^ 

1/2 3.0 xlO^ 

5/2 54.7 xlO~^ 
1/2 974 
5/2 58.3 xlO^^* 
1/2 7.99 xlO^ 
5/2 131 
786 



1 
1 
1 
1 
1 
1 
2 
2 

2 
2 
2 
2 
2 
3/2 



+ 
+ 
+ 

+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



+ 

+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 

+ 
+ 
+ 

+ 



8.186 (10) 
8.287 (10) 
8.436 (10) 

8.470 (10) 



8.680 (7) 

r = 3/2- 

r = 54(6) keV 



8.065 (10) 

r = 3/2+ 

r = 8(3) keV 



8.849 (5) 
r = 10 keV 



8.981 3/2+ 

8.264 5/2" 

8.516 5/2" 

8.16 5/2+ 

8.774 5/2+ 



8.929 5/2+ 



9.099 5/2+ 



9.232 5/2^ 



8.658 9/2" 



no 
a 
ni 



no 
a 
no 
a 
ni 

no 
a 

no 
a 

ni 



no 
a 



ni 



no 
a 



ni 



no 
a 



ni 



no 
a 
ni 

no 



2 
2 
1 
3 
1 
1 
1 
3 
5/2 



5/2 
3/2 
5/2 
1/2 



2 3/2 84.1 
2 5/2 2.18 xlO^ 
2 1/2 5.13 xlO^ 
2 5/2 403 
3/2 19.4 
2 3/2 1.77 xlO"^ 
5/2 1.06 xlO^ 
1/2 881 
5/2 54.4 xlO"^ 
1/2 3.3 xlO^ 
21.3 
1.07 xlO^ 
1.22 xlO^ 
1.86 xlO^ 
16.0 xlO"^ 
2 5/2 0.45 xlO"^ 
2 1/2 23.0 xlO^ 
5/2 217 
2 5/2 32.4 
2 3/2 448 
5/2 10.3 xlO^ 
2 5/2 9.35 
2 1/2 18.7 xlO^ 
5/2 31.5 
2 5/2 768 
2 3/2 1.03 
5/2 4.43 xlO^ 
2 5/2 24.9 xlO"^ 
2 1/2 
5/2 
2 5/2 
2 3/2 
5/2 
2 5/2 
2 1/2 
5/2 
2 5/2 
2 3/2 
5/2 



5/2 
1/2 
5/2 
3/2 
5/2 
1/2 



9.34 xW 

16.8 xlO^ 

1.1 xlO^ 

37.7 

16.7 xlO^ 
816 xlO"-'^ 
55.0 xlO^ 
4.69 xlO^ 
2.92 xlO^ 

20.7 

20.8 xlO^ 
444 xlO"^ 

110 xlO-'^ 

625 xlO"^ 

222 

222 

109 



+ 
+ 
+ 
+ 

+ 
+ 

+ 
+ 
+ 
+ 
+ 

+ 
+ 
+ 
+ 
+ 
+ 
+ 



5.526 (2) [20] 
r = 6 keV 



8.930 (5) 
r = 5 keV 



TABLE III: i^O(a,n)2"Ne and i^O(a,7)^^Ne reaction rates 
(in units of cm^ mol~^ s~^). Above T = 2 GK the recom- 
mended rates are Hauser-Feshbach calculations, below they 
are the sum of the respective low-energy estimates and the 
experimental data. 



{a,n) 
T [GK] NA{av),^p. NA{av)r 



0.1 
0.11 
0.12 
0.13 
0.14 



9.96 xlO 

1.44 xlO 
1.54 XlO 
1.28 XlO 
8.76 XlO 



-20 
-19 
-18 
-18 



5.54 XlO 
6.89 XlO 
6.13 XlO 
4.21 XlO 
2.40 XlO 



-19 
-18 
-17 
-16 



(a, 7) 

Na{cv)exP. NA{crv)re 



9.52 XlO 
4.24 XlO 
4.61 XlO 
1.69 XlO 
2.65 XlO 



-34 
-31 
-28 
-26 



1.67 XlO 
3.37 XlO 
4.17 XlO 
3.50 XlO 
2.16 XlO 



-19 
-18 
-17 
-16 
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{Continued.] 



T [GK] NA{av)e.: 



{a 



0.15 
0.16 
0.17 
0.18 
0.2 
0.25 
0.3 
0.35 
0.4 
0.45 
0.5 
0.6 
0.7 
0.8 
0.9 



5.03 
2.49 
1.09 
4.35 
5.64 
2.66 
4.37 
1.89 
3.28 
3.06 
1.86 
3.06 
2.67 
1.67 
8.53 



xlO"'' 
xlQ-'^ 
xlO"^'^ 
xlO"^^ 
xlO"" 
xl0~" 
xlO"" 
xlO~^ 
xlO"'' 
xlO"^ 
xlO-* 
xlO"^ 
xlO"^ 
xlO"^ 
xlO"^ 



,n) 

NA{ov)r 



N_ 



1.20 
5.44 
2.25 
8.51 
9.33 
1.11 
6.50 
2.14 
3.45 
3.14 
1.89 
3.08 
2.68 
1.67 
8.53 



xlO"''' 
xlO"^'^ 
xlO-i* 
xlO-i* 
xlO"^^ 
X 10-10 
xlO"" 
xlO"^ 
xlO"^ 
xlO"^ 
xlO-* 
xlO"^ 
xlO~^ 
xlO^i 
xlO~i 



(a, 7) 



2.10 
9.58 
2.77 
5.49 
8.69 
7.52 
3.01 
2.10 
4.94 
5.63 
3.88 
6.79 
5.07 
2.24 
7.06 



xlO-"" 
xlO"^^ 
xlO"^i 
xlO-2° 
xlO~i* 
xlO-i'' 
xlO"" 
xlO"^ 
xlO"* 
xlO"^ 
xlO"^ 
xlO"^ 
xlO-* 
xlO^^ 
xlO"^ 



1.04 
4.13 
1.38 
4.03 
2.48 
9.56 
4.35 
1.19 
1.61 
1.30 
7.17 
9.82 
6.51 
2.70 
8.16 



"xicF^ 

xlO'i'^ 
xlO-i^ 
xlO-i'' 
xlO"^^ 
xlO"^^ 
xlO-i° 
xlO"** 
xlO"'' 
xlO"^ 
xlO"^ 
xlO"'' 
xlO--* 
xlO"^ 
xlO"^ 



1.0 
1.25 
1.5 
1.75 
2.0 
2.5 
3.0 
3.5 
4.0 
5.0 
6.0 
7.0 
8.0 
9.0 
10.0 



3.63 
6.65 
5.41 
2.52 
8.10 
4.10 
1.18 
2.47 
4.22 
8.54 
1.30 
1.72 
2.05 
2.32 
2.51 



xlO" 
xlO^ 
xlO^ 
xlO^ 
xlO^ 
xlO* 
xlO^ 
xlO^ 
xlO--' 
xlO^ 
xlO'^ 
xlO*^ 
xlO*^ 
xlO** 
xlO*^ 



xlO" 



3.63 

6.65 

5.41 

2.52 

8.10 

5.01 

1.88 

5.17 

1.14 

3.71 

8.61 xlO 

1.63 

2.71 

4.06 

5.76 



xlO^ 
xlO^ 
xlO^ 
xlO^ 
xlO-* 
xlO^ 
xlO^ 
xlO*' 
xlO*^ 



xlO^ 
xlO'^ 
xlO^ 
xlO'^ 



1.75 xlO"^ 
9.22 xlO~^ 
2.89 xlO"i 
6.77 xlO-i 
1.30 xlO" 
3.34 xlO" 
6.19 xlO" 
9.47 xlO" 

1.27 xlO^ 
1.86 xlO^ 

2.28 xlO^ 
2.56 xlO^ 
2.72 xlO^ 
2.79 xlO^ 
2.81 xlO^ 



1.98 xlO"^ 
9.97 xlO"^ 

3.05 xlO-i 

7.06 xlO-i 
1.30 xlO" 
4.16 xlO" 
9.58 xlO" 
1.82 xlO^ 
3.05 xlO^ 
6.77 xlO^ 
1.22 xlO^ 
1.93 xlO^ 
2.79 xlO^ 
3.79 xlO^ 
4.90 xlO^ 
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